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Abstract 

The electronic states of the two-dimensional Hubbard model are investigated by means of a 4-pole approximation 
within the Composite Operator Method. In addition to the conventional Hubbard operators, we consider other two 
operators, which come from the hierarchy of the equations of motion. These operators carry information regarding 
surrounding spin and charge configurations. By means of this operatorial basis, we can study the physics related to 
the energy scale J = 4t 2 /U. Results present, in addition to the main two bands, a quasi-particle peak at the Fermi 
level, shadow bands and band flatness at (n, 0) point. 
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The main gross feature of the Hubbard model is the 
splitting of the band with the formation of a gap of the 
order U. This physics can be well understood in terms 
of the two Hubbard operators, which mainly describe 
the two Hubbard sub-bands. However, in order to catch 
low-energy features related to the energy scale J = 
4t 2 /[/, it is necessary to take into account higher-order 
operators, which carry information regarding nearest- 
neighbor spin and charge correlation effects. Along this 
line, we analyze the Hubbard model by means of a 
4-pole approximation within the Composite Operator 
Method [1,2,3,4]. In addition to the conventional Hub- 
bard operators, we consider other two operators that 
come from the hierarchy of the equations of motion. 
The resulting band structure will be much richer than 
the 2-pole one. 

The two-dimensional Hubbard Hamiltonian reads as 
follows, 

h = X)(*u - /^u)4«<v(j) + 5> T (i)»u(0, (1) 
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where c£(i) and Co-(i) are creation and annihilation op- 
erators of electrons with spin a at the site i [i = (i, t)], 
respectively. n CT (i) = cJ.(i)c CT (i), [X is the chemical po- 
tential, tij = — 4tay, a[k] = ^[ay] = ^(cos(k x a) + 
cos(fc !/ a)), a is the lattice constant, T is the Fourier 
transform, U is the on-site Coulomb repulsion. We de- 
fine the following operatorial basis, 



(2) 



£o-(i) = Co-(i) (1 — n_ CT (i)) and r] a (i) = c CT (i)ra_ CT (i) de- 
scribe the transitions n(i) = <-> 1 and 1 <-»> 2, respec- 
tively. The equations of motion of ipA<r(i) give 



g 
d 

i-^Vcr(i) = (-A* + U)ri<j{i) + 4t7r CT (i). 
at 



(3) 



where 7r CT (z) = — n_ CT (i)c"(i) + cL CT (i)c CT (i) c -<r(i) + 
c CT (i)c" t CT (i)c_ CT (i) with c"(i) = V\ a y c CT (j,t). Now, we 
divide 7r CT (i) into two operators 7T CT (i) = £sa(i) + rj aa {i) 
in the same manner as we have done with c„(i) = 
£<r(i) + f)<r(i)- Then, we define a new operator set 
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with 



ri,a(i) 



(4) 



£ SCT (i) = -n-„{i)£(i) + cl CT (i)c CT (i)£™ CT (i) 
+c CT (i)j7" t CT (i)c_ CT (i) 

Vsa{i) = -n- a (i)r)"(i) + c^_, y {i)c CT {i)7]" a {i) 



(5) 



It is worth mentioning that ipBa(i) describes two-site 
composite excitations [5] as they carry information of 
surrounding spin and charge configurations, and are 
eigenoperators of the interaction term as £ CT (*) and 

In the present paper, we choose as operatorial basis 



■4>Acr{i) 
1pBcr(i) 



(6) 



Within the Composite Operator Method, once we 
choose a n-component spinorial basis ip, the equations 
of motion take the general form 



d 



(7) 



where e(k) = m(k)J 1 (k) is n x n matrix as J(k) = 

t), (j, *)}> and m(k) = T({i£ t ^(i, t), (j, *)}>. 
If we neglect Sj we obtain a pole structure for the 
Green's function 



(8) 



The 2-pole solution with ip — ipA have been discussed 
in detail in Refs. [2,3]. In the present paper, we analyze 
a 4-pole solution by means of the basis (6). According 
to this, we have 



\l A B(k) /Bfl(k) J 



(9) 



hermiticy. For e_gs(k), we use the simplified equations 
of motion discussed in Ref. [5], 



i-p^r)sa{i) ^ {-fJ, + U)r] sa {i) + At 



(11) 



where we neglect terms that give diffusion processes 
over three sites. This procedure simply gives ess (k) by 
inspection. The correlation functions appearing in the 
I and e are self-consistently determined by means of 
Green's function and the local algebra constraints [4]. 

In Fig. 1, we present the density of states and the 
dispersion relation. It is worth noticing that the above 
formulation is applicable in any dimension. In addition 
to the main Hubbard band structure, our results show 
a coherent peak around the Fermi level and shadow 
bands originated by the antiferromagnetic correlations. 
The details of formulae and a more extended analysis 
will be presented elsewhere. 
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e(k) = 



€AA (k) CAs(k) 

esA(k) ess(k) 



(10) 



The main difficulty of the present formulation is the 
evaluation of Jss(k) and eBs(k). The other elements 
can be simply computed either explicitly or in terms 
of those latter. For Jss(k), we calculate anticommuta- 
tors explicitly and decouple the higher-order correla- 
tion functions preserving particle-hole symmetry and 
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Fig. 1. The density of states (top) and the dispersion relation 
(bottom) for U = 8i, n = 0.9 and T = O.Olt. The width of 
dispersion line gives the intensity of peak. 
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